Analysis of surface plasmon excitation at terahertz frequencies 
with highly-doped graphene sheets via attenuated total reflection 

Choon How Gan 
Electronics and Photonics Department, 
A*STAR Institute of High Performance Computing, 138632, Singapore 

Abstract 

Excitation of surface plasmons supported by doped graphene sheets at terahertz frequencies 
is investigated numerically. To alleviate the momentum mismatch between the highly-confined 
plasmon modes and the incident radiation, it is proposed to increase the surface conductivity of 
graphene through high doping levels or with few-layer graphene. For currently achievable doping 
levels, our analysis shows that surface plasmons on monolayer graphene may be excited at oper- 
ating frequencies up to about 10 THz (~ 41.3 meV) with a high-index coupling prism, and higher 
frequencies/energies are possible for few-layer graphene. These highly-confined surface modes are 
promising for sensing and waveguiding applications in the terahertz regime. 
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Electromagnetic fields associated with surface plasmons (SPs) guided by highly- 
conductive metal surfaces at midinfrared (MIR) and terahertz (THz) frequencies are ex- 
tremely unconfined; they are plane-wave like and spread ubiquitously in the interfacing 
dielectric medium. This poorly-confined nature of SPs propagating along near-perfect con- 
ductors may be understood by inspection of its wavenumber k sp . At a metal-dielectric 
interface^, k sp = k €de m / (e^ + e m ) , where k is the wavenumber in free-space, and 
and e m are the dielectric constants of the dielectric and metal, respectively. For the sim- 
ple half-space problem, the transverse wavenumber of the SP in the dielectric is given by 
^zd = ^o €d ~ k 2 sp = kle^/(ed + e m ). At visible and near-infrared frequencies, |9fte m | >> Qe m 
(typical of the noble metals), the imaginary part of k z d dominates and the field is well- 
confined. At THz frequencies however, |3?e m | and both large with comparable 
magnitudes. Typically |e m | >> e^, k sp ~ koy/ea, and k 2 zd ~ /cQe^/e m , the imaginary part of 
k Z( i is relatively small leading to a long decay length in the dielectric. 

To circumvent the poor confinement at long wavelengths, tunable SP-like modes sup- 
ported on metal surfaces corrugated with nanoholes or nanorods, and on multi-layered meta- 
material systems have been proposed^-. More recently, studies have shown that graphene, 
a single layer of carbon atoms gathered in honeycomb lattice, can also support well-confined 
SP modes at MIR and THz frequencies^-. As its electronic transport properties can be read- 
ily tuned by the application of a gate voltage*^, graphene structures present an attractive 
alternative as platforms for supporting SPs. However, due to the relatively large momentum 
mismatch between the SPs and the light photons, their excitation on graphene remains a 
challenging issue. Proposed approaches to overcome the mismatch include the patterning 
of graphene into an array of nanoribbons^i, and probing with a nano-sized metalized tip in 
near-field nanoscopy^ 1 ^. In this letter, we investigate the excitation of SPs on highly-doped 
graphene sheets with attenuated total reflection (ATR) via the Otto geometryi^. As shown 
in Fig. 1, graphene is deposited on a substrate of refractive index m, which is separated 
from the coupling prism (n p ) by a small gap d {n-z). In the following, the substrate is taken 
to be nonpolar so that effects of remote phonon scattering may be neglected^. Unlike in 
metallic thin films where the SP dispersion splits into two branches, monolayer graphene 
(MLG) sandwiched between two dielectric media supports a single bound SP mode^. For 
well-confined SP modes supported on MLG, the associated magnetic field (H) and normal 
component of the electric field (E z ) change sign across the graphene sheet on account of the 
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FIG. 1: Exciting SPs on MLG/FLG with the Otto geometry. 

continuity of the tangential electric field in the two neighboring media with positive permit- 
tivity values. In addition to highly-doped MLG, few-layer graphene (FLG) will be considered 
in our attempt to reduce the momentum mismatch. By taking individual graphene sheet as 
a non-interacting monolayer, the optical conductivity of the FLG^>^ is Na, where N is the 
number of layers (N < 6) and a is the conductivity of the MLG. As will be shown shortly, 
increasing the conductivity effectively reduces the SP wavenumber. 

Let us model each graphene monolayer as a surface conducting sheet. Near to the THz 
regime, intraband scattering dominates in highly-doped graphene and its conductivity takes 
on a Drude-like^- 1 ^ form a ~ ie 2 fj,/[7rh 2 (uj + iT -1 )]. Here \x ~ fivp^fW^n is the Fermi-level, 
vp ~ 10 6 m/s is the Fermi velocity, N c is the carrier density, u = k c is the angular frequency 
with c the speed of light in vacuum, r is a phenomenological electron relaxation time, and 
the local limit u » r~ l ,k sp Vf is assumed. By matching the boundary conditions for the 
?ii-MLG-n2 system (Fig. 1), the SP dispersion may be derived as^ 1 ^ 
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If §lk sp » ^Sk sp and \k sp \ 2 » kf, where fc, = njfc and % = {1, 2}, the transverse wavenum- 
ber of the SP is k zi = U/k 2 — kf w ik sp . Substituting k zi = ik sp into the above dispersion 



relation yields the approximation^^ 



(2) 



From this expression, it is readily seen that the SP wavenumber decreases with increasing \x 
or a. Figure 2 shows the SP dispersion for the MLG as a function of frequency / in the THz 
regime and as a function of the Fermi-Level \x. For all following calculations, r is taken to 



be 1 ps to account for scattering loss from acoustic phonons^ 1 ^ -, and the refractive indices 
rix and ri2 are taken to be 1.6 and 1, respectively. In Fig. 2(a), the Fermi-level is taken to be 
/i = 0.9 eV, which corresponds to iV c « 6 x 10 13 cm" 2 , within the doping levels achievable 
in experiments 2 ^ 2 -. The decrease in the effective mode index of the SP (n sp = k sp /ko) 
as the frequency decreases is evident. The data also shows that SPs can propagate for 
hundreds of microns (L = A/ (47r9 ; n S p) ) in the THz regime with effective wavelengths 2-6 
times smaller than the free-space wavelength A = 2tt jk§. For / = 5 THz, the dependence of 
the SP dispersion on \x is shown in Fig. 2(b), from which the decrease of n sp with increased 
fi is clearly observed. The relations §tn sp oc u,fi~ l and L sp oc derived from the 

approximate form Eq. ([2]) (dashed curves in Fig. 2) is in good agreement with the general 
trends of the SP dispersion. However, as seen from Fig. 2, the approximated dispersion 
relation (J2J) starts to deviate significantly from the exact form (Eq. (0Q)) for frequencies 
/ < 3 THz or Fermi- levels fi > 1.6 eV, especially for the imaginary part of the dispersion. 
Beyond these limits, the imaginary part of k zi is no longer the dominant term, leading to 
the inaccuracies in the approximation ([2]). 

The excitation of SPs with the MLG or FLG in the Otto geometry is demonstrated 
with the calculated angular reflection spectra in Fig. 3, taking the operating frequency to 
be / = 5 THz (hu = 20.7 meV). In our calculations performed with the transfer-matrix 
method^ 3 -, it is taken that each graphene sheet has effective thickness^ t g ~ 0.5 nm and an 
equivalent dielectric constant e g = 1 + i<7/W t 9 . For highly-doped MLG, \3te g \ » Se 5 in 
the local limit (u » r" 1 ) since Scr/3ftcr « ujt, providing an intuitive insight to why highly- 
confined SP modes may be supported by graphene even at THz fequencies. The variation of 
the SP dispersion with the number of graphene sheets is shown in Fig. 3(a), with the Fermi- 
level /i taken to be 0.9 eV. By increasing N gradually from 1 to 5, the SP propagation 
distance increases from about 60 /xm to more than 1 mm while its effective wavelength 
increases from roughly A/3.2 to A/1.6. Figure 3(b) shows that a significant portion of the 
SP field extends to the neighboring media. Of particular interest is the penetration into 
the gap region (71,2) as both the width and the angular position of the resonance dip is 
modified in the presence of the coupling prisrn^. The width of the resonance, which is 
roughly proportional to ^sk sp , is broadened by radiative damping of SPs that couple back 
to the prism as freely-propagating photons. For narrower gaps d or where the SP fields 
are more intense, stronger coupling with the prism also leads to a more profound angular 
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FIG. 2: SP dispersion on MLG as a function of (a) frequency, and (b) Fermi- level. Solid curves 
are results from Eq. (pQ), and dashed curves are calculated from the approximation in Eq. ([2]). 

shift of the ATR minimum from the angle sm~ {n sp /n p ). It is observed from Fig. 3(b) 
that for gap distances d < 10 /im, significant coupling may occur between the SP and the 
prism. Calculations of the transverse decay length of the evanescent field at the prism/gap 
interface show that the gap distance should be limited to d < 5 /im so that it may couple 
to the surface mode. Therefore, it may be expected that the radiative damping of the SP 
can affect the ATR minimum non-negligibly. 

The systematic shift in the ATR minimum from </? = 55.8° for iV = 1 to ip = 24.8° for 
iV = 5 is shown in Fig. 3(c), taking the refractive index n p ~ 4.0 for a Germanium^ prism. 
The resonance dip is clearly resolved in each case, with a difference of 0.8° between N = 4 
and N = 5. For the MLG (N = 1), the reflection minimum reaches 50%, and the width of 
the resonance Aip = 4.52°. For N > 2, the strength of the SP field that extends into the gap 
decreases, and sharp narrow dips are obtained with the reduced radiation damping of the 
SP mode. The shift of the ATR minimum from the angle sin -1 (n sp /n p ) is 3.45° for the case 
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N = 1, and gradually decreases for increasing N to 0.13° for N = 5. Similar trends in the 
strength and width of the resonance dips are observed in the reflection spectra of Fig. 3(d) 
for MLG with its Fermi-level increased from \x = 0.8 eV to 1.6 eV (N c m 2 x 10 14 cm -2 ). 
From the simulated reflection spectra, it is observed that the narrowest and most pronounced 
resonance dips occur for the least lossy SP modes, i.e. for increasing number of layers iV or 
Fermi- level [i. 

It is reasonable to pose the question on the range of frequencies and doping levels for 
which SPs may be excited on graphene with the Otto configuration. Let us revert to the 
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FIG. 3: SP dispersion and simulated angular reflection spectra with MLG and FLG for / = 
5 THz (fruj = 20.7 meV). (a) SPP dispersion for N = 1 to 5. (b) Comparison of the normal electric 
field component E z for N = 1 and N = 5. (c) Spectra for N = 1 (top) to N = 5 (bottom). In (a) 
to (c), the Fermi-Level \i is taken to be 0.9 eV. (d) Spectra for \i = {0.8, 1.0, 1.2, 1.4, 1.6} eV (top 
to bottom) for MLG (N = 1). For all cases in (c) and (d), the gap distance d was varied between 
3 — 6 /j,m. 
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approximate dispersion relation (TJJ) to provide an estimate. Imposing the constraint $tn sp < 
n p leads to the condition 

A > 8{n\ + n 2 2 )/n p , (3) 

where 5 = 2eo(vr/ic) 2 /(e 2 /x). Let us take that the refractive index of the substrate is greater 
than unity, i.e., n\ > 1. The limit of the inequality expressed in condition 03]) is plotted for 
different values of Tip clS db function of the Fermi-level in Fig. 4, taking the gap region to be 
free-space (n 2 = 1). The domain above each curve is where the inequality is satisfied. For 
current achievable levels of doping for graphene 2 ^ 12 ^, Fig. 4 shows that the highest frequency 
for which SPs on MLG may be excited through ATR with high-index coupling prisms is 
~ 10 THz (A ~ 30 yum). For excitation at longer wavelengths (provided A << 2ir cr), heavy 
doping to increase the Fermi-level becomes essential. The condition ([3]) in terms of photon 
energy (in units of eV) is hu < efm p / [eotr hc{n\ + n|)]. As the inequality holds for MLG, it 
suggests that excitation frequencies / > 10 THz (Tvjj > 41.3 meV) may also be considered 
in the case of FLG. 

In conclusion, we have demonstrated numerically the excitation of SPs at THz frequen- 
cies with highly-doped MLG and FLG via the Otto geometry. In contrast to the noble 
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FIG. 4: Estimating the spectral range for which SPs may be excited on MLG with the inequality 
A > 5{n\ + n\)/n p (Eq. ([3])), taking the gap region to be free-space (ri2 = 1) and n\ > 1. The 
inequality is satisfied in the domain above each curve. The photon energy (in units of meV) 
corresponding to the wavelength is shown on the top horizontal axis. 
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metals, the SP modes supported on graphene sheets remain highly-confined even in the THz 
regime. This strong confinement may be explained by considering an equivalent dielectric 
constant for graphene, from which it is apparent that \dte g \ » Qeg. With rapid advances 
in the fabrication of MLG and FLG, it is foreseeable that the above simulated results may 
be realized in a future experiment. For FLG, in case each of the graphene sheet is separated 
by a transparent buffer layer (~ 20 nm thick) as in Ref. [17], our calculations show that the 
effect on the position, width, and magnitude on the resonance dips is only marginal. As an 
alternative, the Kretschmann configuration may be employed if the substrate containing the 
graphene sample may be directly deposited onto the coupling prism. If high-index material 
at THz^i becomes available, higher-index prisms can be employed to excite SPs on MLG 
for frequencies beyond 10 THz. High-index prisms are also useful to shift the resonance to 
smaller angles where the transverse decay of the evanescent fields are longer, allowing for 
larger gap distances d. Apart from surface plasmon resonance sensors, the proposed con- 
figuration may be employed for determining the thickness of graphene samples containing 
MLG and FLG^>2& Finally, it is worth noting that in contrast to the H— symmetric mode 
supported in thin metal films, the field of the H— antisymmetric SP mode supported on 
graphene decays more quickly in the medium of lower refractive index (see Fig. 3(b)). As 
such, the SP mode may be considered for improving the confinement factor in THz waveg- 
uides, for instance in semi-insulating surface plasmon waveguides of the quantum cascade 
laser—. 
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